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1 Introduction 



Precise experimental measurements in high energy physics require, in its theoretical counterpart, the 
development of techniques for the evaluation of analytic objects associated with the corresponding 
Feynman diagrams. These techniques have lately emphasized the automatization of calculations of 
multiscale, multiloop diagrams. 

Modern numerical methods for the evaluation of Feynman diagrams benefit from analytical 
techniques employed as preliminary work to detect the presence of divergences. Recent advances 
include a method based on the Bernstein-Tkachov theorem for the corrections of one and two loop 
diagrams p], a method based on sector-decomposition, as applied in [2] to a variety of diagrams 
and to the evaluation of integrals in phase space, and a third example is the one developed in [3], 
which contains numerical evaluations of two-point functions by differential equation methods. 

New analytic methods include techniques to reduce Feynman diagrams to a small number of 
scalar integrals, such as integration by parts [4], the use of Lorentz invariance [5] and the use 
of symmetry 6 . Methods for the evaluation of scalar integrals include expansion by regions [TJ, 
Mellin-Barnes transforms |5], relations among integrals in different dimensions [5] and differential 
equations (It)] . 

This paper contains examples of an alternative method for the evaluation of some Feynman 
diagrams. The ideas are based on the classical Ramanujan Master Theorem (RMT), which is a 
favorite technique employed by the well-known mathematician to evaluate definite integrals. The 
theoretical aspects of this method are presented in |12j and a collection of modern examples is given 
in [13]. The application of this technique has been illustrated in [IT] with the evaluation of some 
multidimensional integrals obtained by the Schwinger parametrization of Feynman diagrams. The 
goal of the present work is to use this technique to evaluate integrals associated to two and three 
loop diagrams. The method works for a large variety of definite integrals and our first example 
illustrates this by computing the Mellin transform of a Bessel function. The next example illustrates 
the evaluation a multidimensional integral corresponding to the massless bubble Feynman diagram. 
Then we provide a generalization of RMT to multiple integrals, required for multiloop calculations. 
Further applications will be described in future work. 



2 Ramanujan's Master Theorem (RMT) and its generaliza- 
tion : The formalism 

Integrals of the form J °° dx x v ~ x f(x) may be evaluated by one of Ramanujan's favorite tools; the 
so-called Ramanujan Master Theorem. It states that if f{x) admits a series expansion of the form 

oo ( \n 

' 71! 

in a neighborhood of x — 0, with /(0) = <p(0) 7^ 0, then 

00 

J dxx"- 1 f{x)=T{u) V {-v). (2) 


The integral is the Mellin transform of f(x) and the term <p{—v) requires an extension of the 
function (p, initially defined only for v G N. Details on the natural unique extension of ip are given 
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in [T3]. Observe that, for v > 0, the condition <p(0) ^ guarantees the convergence of the integral 
near x — 0. The proof of Ramanujan Master Theorem and the precise conditions for its application 
appear in Hardy |12j . The reader will find in |13j many other examples. 



2.1 The Mellin transform of a Bessel function 

The first illustration of RMT involves an integral containing the classical Bessel function, given in 
its hypergeometric form by 



1 



2 ) r(l + a) 



1 + a 



Here qFi is the hypergeometric function defined by 



E 



1 x r ' 



(a)„ n 



using the Pochhammer symbol 



{a) n 



T(a + n) 



The integral to be evaluated here is 



oo 

I = J dx X^Jo, (s/x) , 



and can be expressed as 



I = I dx x?- 1 ( v± 



dx ~ ( - ir 

X n =0 



1 



E 



(-1)" 1 x n 



r (1 + a) „tl o n! (1 + a). 4™ 



2 Q + 2n r(l + a + n) 



The data matches the statement of RMT with 

(p(n) 



Therefore, the evaluation 



2«+ 2 »r(l + a + n)' 

r(-n») 

2"+ 2 «*r(l + a + n*) 



(3) 



(4) 



(5) 



(6) 



(7) 



(8) 



(9) 



is obtained directly from RMT. Here n* = — (8 + ^ ). For information to be used in the latter 
sections, observe that n* is the solution to 



n + (3+ f = 0. 



(10) 
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It follows that 



oo 

/ 



dx x^J a = 2^ P^fi v (11) 



This result appears as entry 6.561.14 in the table of integrals |15) . 

2.2 A second example : the Feynman diagram of a bubble 




(12) 



The next example illustrates the evaluation a multidimensional integral corresponding to the mass- 
less bubble Feynman diagram shown in (fT2")l . The result is well-known [14] . In momentum space 
the corresponding integral is given by 



G = 



d D q 



■ I 1 ■> r nl Q 2 ' (13) 

k 2 r \(p-if 



where the parameters {a,} are arbitrary. The Schwinger representation corresponding to this dia- 
gram produces 

f-ir* 77 exp (-7T^ p2 

G= j r , ; / / 'A-/// -r" '//"• : V +V D 7 - (14) 

r(ai)r(a 2 ) J 7 (x + y)^" 

In order to apply RMT to evaluate this integral, each part of the integrand is expanded in a 
Taylor series. In situations where options are available, the optimal course of action seems to be 
to minimize the number of expansions. At the moment this is at the heuristic level. In the present 
example, it is convenient to first expand the exponential function 

V x + v J i x + y) 

to produce the expression 



D OO OO 



G = < / fdxdy x^-Y*- 1 E <=jf (PT f % +n - (16) 

1 (ax) 1 (a,2) J J „=o (x + y) 2 + 





The final step is to expand the binomial term (x + y)~ D / 2 ~ n in the form 



oo 

£<#(# + n) k x-^— V, (17) 



D 

(x + y) 2 +n k=Q 



where the Pochhammer term is 
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(#+»), 



*) 



r(f + „) 



Replacing in (JTSJ) gives 



oo oo 



]j oo oo 

(V = J^G ; I I — — E E +»). 



r(ai)r(a 2 ) J J x y k=Qn =o 
o o 

and x — ► 1 /x produces the alternative expression 



G 



U OO OO 

(-l) -ir /" f dxdy 



r(ai)r(a 2 ) J J x y k=Qn =o 
o o 



V V (-1)" (-!)*■ 



There are several options to employ RMT to evaluate this integral : 

a) Evaluate the integral in the x- variable and the expansion in the index k 



oo 

/ 



dx... E • 

fc=0 



{-X? 

' k\ 



b) Evaluate the integral in the y-variable jointly with the series with index k 



dy... E ■ 

fc=0 



-y) k 



k\ 



c) Evaluate the integral in the y- variable jointly with the series with index n 



dy... E 



i-vY 



i=0 n - 


It is now shown that each of these options produces the same result. 



2.3 Solution with option (a) 

The hrst procedure to evaluate G is given by 



D OO 

(-1) T f d y 



T(a 1 )T(a 2 ) J y 



where <p(k) is dehned by 



TV. 

n=0 



The value of (|2Tj) is given by RMT as 



D OO 



G[a) = J^- } J f r(-k*)v (k*) , with k* = -f + fll 



o 



— -r( T - oi) / — E HtMp ) {T + n ) ai -% v 





r(oi)r(a 2 ) 



ai+o 2 



r(ai)r(a 2 ) y y »=o r(^ + n) 



o 



RMT is applied one more time to the last integral to obtain 

■* = r (ai )r(a 2 ) r( ^ ^ ai)r( - n } ( p ) r(f + n-) ' 

with n* = — oi — a 2 + y. Therefore, option (a) gives the value of G as 

^ , n _ d , 2 °- ai -a 2 r( ai +a 2 - f )r(f -ai)r(f -a 2 ) 

G (a) = 2 J 



r(ai)r(a 2 )r (D - oi - a 2 ) 



2.4 Solution with option (b) 

A similar argument now yields 



r(oi)r(a 2 ) J x J y k=0 h - 
o L o 



o 

with 

oo 

^) = £^(pT(f+«)* 

Therefore 



fe— ai+- 



G 



G ib) = f gf^y/^r(- fc >(^ ) , with fc * = -n- 

o 



a 2 







— n — a 2 



j,— n — a 2 — 01 + 



(-i) f Jdx^i-ir r ,vn r(n + a 2 )r(f - ffl2 ) 
r( ai )r(a 2 )i x ^ - ! lPj r(f + n) 

The change of variables a; 1 /x gives 



x — n-a 2 -ai + -f 



r (-1) 

G (6) = — - 



r( ai )r(a 2 )i x ^ m r(f + n) 



and RMT produces 



with n* = — a 2 — 



r - (-!)' f n „*) (n 2^ r(n* + a 2 )r(f -a 2 ) 

-a 2 — ai + -j. Therefore, the value of G obtained from option (b) is 

G{b) = {p ] r( fll )r(a J )r(fl- aj - ai ) 

2.5 Solution with option (c) 

Proceeding as before produces 

D OO 

(-1) ^ /■ 
G (c) - — — 



with 



r(a 1 )r(a 2 )7 x i y io " ! ^ ( j y 

o L o 

OO 



This yields 



with n* = —k — a 2 . Then 



n 



D °° 



o 

with the equivalent form 



r 2 rfg n\( dx ^ r^\- fc -^ r(fc + a 2 ) „ fc _ ai+ g , . 

G ^ = r(a 1 )r(a 2 ) r (-- a2 ) J /- fe ? ^( p ) r(#-fc-a 2 )" " (36) 

o 

An application of RMT gives 

Mr* t,/d w 2\-fc*-«2 r(-fc*)r(r + a 2 ) _ £> 
G ^> = r (ai) r M r (t - «*) (p 2 ) r(#- fc *- a2 ) Wlth fc = fll " * ■ (37) 

or equivalently 

r ( 1 V-3 ^ r (f - Q 2) F (f - a F ( Q 1 + fl 2 - f ) /ooN 

<^( C ) = (-1) (P J ^ — r-p7 — vFTn ^ ' ^> 

T (oi) r (o 2 ) F (£> — ai — a 2 ) 

Observe that each option produces the same value for the integral G; that is, G( a ) = G^ = G( c ) . 

3 Generalization of RMT to multiple integrals. Multiloop 
calculations 

This section discusses a generalization of RMT to multidimensional integrals of the form 

oo oo 

/= [dxix?- 1 ... f dx N x'X' 1 f(x u ...,x N ). (39) 



As in the one-dimensional case, the function / = / (xi, Xn) is assumed to admit a Taylor 
expansion given by 



oo oo 



i 1= o i N =o (40) 

v ~.Ollil + -"+OlJvijV+6l ^INlh+.-.+aNNlN+bN 
A ... , 

so that I is expressed as 

oo oo 

f dxx f dx N °° °° (_i)i N 
J x\ j xn i n= q 



(41) 



v ™anii+...+aiAr/jv+bi ^aivi^i + .-.-f-aivivJiv+friv 
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with bi = v% + bi (i = 1, N). 

Applying RMT in iterative manner produces 



/ = 



1 



|det (A) 

where det (A) is the determinant 



an 



aiN 



det (A) = 

O-Nl ■ ■ ■ O-NN 

and the variables I* (i = 1, ...,N) are solutions of the linear system 

anh + ... + duvijv + b\ =0 

(iNili + ••• + unnIn + bpf = 0. 



(42) 



(43) 



(44) 



Details of the proof of this result appear in [T3]. This procedure will be called the Generalized 
Ramanujan's Master Theorem (GRMT). 

4 Applications 

4.1 Massive sunset diagram 

The integral evaluated first is associated to the diagram shown in the figure. In momentum space, 
the integral is given by 



G = 



d u qi d u q 2 1 



1 



i7r D/2 ijr D/2 [ g 2 _ M 2y 



[(P + 92 



|21 a 3 



(45) 



G 

















(46) 



In terms of the Schwinger parametrization, G becomes 

exp (xiM 2 ) exp 



G = 



(-1) 



—D 



OO OO OO 



r(oi)r(a 2 )r(a 3 ) 



x x x 2 x 3 



X\X 2 + X1X3 + X2X3 



n n n {xix 2 + X1X3 + x 2 x 3 y 



(47) 



where dlzf = dx\ dx 2 dx% . In order to illustrate the power of the GRMT method, 

the special physical case p 2 = M 2 is considered. The integral is now written as 
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G = 



(-1) 



-D 



oo oo oo 



T(a 1 )T{a 2 )T{a 3 ) 



exp 



x\ (x 2 + x 3 ) 

X\ (x 2 + X 3 ) + X2X3 



000 

The expansion of the exponential function gives 



exp 



o\ (x 2 + X3) 



^Xi (x 2 + x 3 ) + x 2 x 3 
and this produces 



= E ^ (~ M2 ) 

m=0 



[xi (x 2 + x 3 ) + x 2 x 3 ] 2 

2ni / 1 \«i 
ni X l [X 2 + X 3 ) 



[Xl (X2 + X 3 ) + X 2 X 3 ] 



(48) 



(49) 



G 



(-1) 



-D 



00 00 00 



V (ai)Y (a 2 )T (a 3 ) j j j ni=0 
000 



[xi (x 2 + x 3 ) + x 2 x 3 ]~ +ni 



(50) 



The binomial theorem is employed next to expand the integrand. This can be done via the hyper- 
geometric representation 



or by the explicit formula 



The result is 



{l + x) a = !F 



n=0 ^ ' 



[xi (ar 2 + x 3 ) + 332X3]' 



r+ni 



E (-i)" 2 r (f + n i 
n 2 ! ~T ( D 



«2=0 



r(f + nij 



X X x (X2 + x 3 j 



° 2 ^3 1 



and (I57J)) is now written as 



G 



(-1) 



-D 



00 00 00 



^ 00 00 

r(ai)r(a 2 )r(a 3 ) 7 7 7 n\— ri2— 





ni! 712! 



* r(#+nr) 

Only the binomial (x 2 + x 3 )~%~ n2 needs to be expanded. This is done as before to produce 



(51) 



(52) 



(53) 



(54) 



(x 2 + x 3 y 



E 

"3=0 



(_!)», T(%+n 2 + n 3 ) 



r(# + n 2 ) 



--"2 ""3^3 



(55) 
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The point has been reached for a direct application of GRMT to evaluate 

OO CO oo 

( -, III dxl dx2 dx3 

Xi x 2 x 3 



oo oo oo 



x E E E ^ 

ni— ri2— ri3— 



(-1)" 2 (-1)" 3 ^ a i~f +ni-n 2 Q3— r -"3„a 3 +n 2 +n 3 



n 2 ! «3 ! 1 



(-1) ° T(f +n 1+ n 2 )r(f +n 2 +n 3 ) , ni 

r( ai )r(a 2 )r( fl3 ) r(f + m) r(f + n 2 ) 1 ' 

It follows that the integral G is given by 



(56) 



G = 



(-1) 



r(oi)r(a 2 )r(o 3 ) 



T(-nl)T(-n* 2 )T(-nt) 



r(f + n 



i-r"5)r(f 



+ n*) 



r(f + „;) 



r(f + n|) 



(-M 2 )" 1 , 



(57) 

where the values assigned to the indices {n*} are the unique solution to the linear system obtained 
from (|56|) . In detail, 



m - n 2 = —a\ 
n 3 — a 2 - y , 
n 2 + n 3 = -a 3 , 



D 
2 ' 



with solution 



Replacing in (f5Tf yields 

g= (-iy 







= -a\ 


- a% — CL3 




< 




= —a-2 


-a 3 + #, 








_ D 
2 


+ a 2 . 




D r («l 


+ a 2 


+ a 3 - 


D) r (a 2 4 


a 3 -f)r(f-a 2 ) 



r( ai )r(a 2 )r(a 3 ) 



r (f - a 3 ) r (2£> - ai - 2a 2 - 2a 3 ) 
ai - a 2 - a 3 ) r (D - a 2 - a 3 ) 



r( 



3D 



2 \ D — ai — 0.2—0.3 



4.2 Massless three loops ladder diagram 



(58) 



(59) 



(60) 



The next example gives the evaluation of the integral associated to the diagram seen in Eq. (|61j) . 
To illustrate the method in a relatively simple situation, the conditions Pf = for 1 < i < 4 and 
,s = arc imposed. 
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G 




(61) 



1' 



SP 3 



The parametric representation of this diagram is given by the integral 



G = 



-l)- v /■ ^ cxp ( — — *) 



r(ai) 



...T(a 10 )J J 







(7 2 



(62) 



where eTaf = Ilj=i x j 3 1 ana ^ tne polynomial t/, written in a form adapted to the application 
of GRMT, is given by 



U = x 5 (x 7 + fi) (f 2 + x 4 ) + x e (x 7 + fi) (f 2 + x 4 ) + x 4 (x 7 + fi) f 2 + x 7 (f 2 + x 4 ) fi, (63) 



with 



fl = X 8 + Xg + X 10 , 



h = xi + x 2 + x 3 
Expanding the exponential term produces 

3D CO OO 



G = 



(-ir 



r(ai) ...T(a W ) J J ni=0 nil 





n i n\ m n\ ^ni 

(ty 



JJ-+ni 



and the polynomial U is expanded using the multinomial theorem 
Od + .-. + Zk-i+a*) = 2^-2^ ^TT 



' ni!...nfc_i!r(l + a-ni - ... - n k -i) 



v ~"i »-i «-m-...-nn 
x l ••• x fe-l x fe 



written in a form adapted to GRMT 



(64) 



(65) 



(66) 



OO OO 



m=0 n fc _i=0 



r (-a + m + ... + n k -i) 



r(-a) 



(67) 



v ~ni o-m-...-n fc -i 
x l ••• x fc-l x fe 
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These expansions produce 



Similarly, 



and 



oo oo oo -n / -D i i i i \ 

U-^-m = V V V (-1)" 2 (-1)" 3 (-1)" 4 r (t + n i + n 2 + "3 + n 4 ) 

' - ' ' ' ' n 2 \ n 3 \ n 4 ! V(—-i-rt-,\ 
n 2 =0n 3 =0n 4 =0 V 2 111 ) 

p-§--ni-n2-n3-n4 fn2 , f ,n 2 +n 3 +n 4 , f \-%--ni-n 2 

X tj t 2 {X 7 + ti) {X4 + T2) 

v „n 2 na„m ~§ -n 1 -n 2 -n 3 -n 4 

y\ Jb ^ 5 Q 7 * 



f-r 1 f \"2+n 3 +n4 _ (-1)" 5 r (~ n 2 - n 3 - n 4 + n 5 ) „ 2+ „ 3+ „ 4 _ ra5 



r-r -u f _ (-i)"» r (f + ni + n 2 + n 6 ) g _ ni _ n2 _ n6 

( 4 ^ " r(# + n 1+ n 2 ) * 4 f2 

The result is 

3D OO 

G r( ai )...T(a W )J -J d nW'nho " l! ^ " 6! W 


r (-j + ni + n 2 + ra 3 + ra 4 ) r (-ra 2 - n 3 - n 4 + n 5 ) r + ni + n 2 + n 6 ) 

r(f + m) T(-n 2 -n 3 -n 4 ) r(f + ni +n 2 ) 

v T H „tii + "2 +n 6 _n 3 _n 4 f _n 2 ""a - n 4+"s ni „- -J -ni -n 5f - # -«l -"6 

A X-^ 10 1 2 

Finally, the powers of fi and f 2 are expanded in the form 

-l)" 7 (-1)" 8 r (f + m + n 5 + n 7 + n 8 ) 



00 00 



ffT -n 1 -n. = 



of^o n ? ! ng! r(f+n 1+ n 5 ) 



n 7 ^.n 8 ---^1 -n 5 -n 7 -n 8 



S\ *Vo <Xj Q X' 



9 -no 



and 



^-f-ni-ne _ (-1)"" (-l)"m 1 j~ + n l + n 6 + "9 + "-10 J 



ng— nio— 



" 9! " lo! r(f+n!+n 6 ) 



v ."s»»io-"T -ni-ne-ng-nio 

At this point the expression for G is in the form required to apply GRMT : 
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3D OO 

(-1) ; jdxi [ten g g (^i...(^^ (ni ,..., nio) 



G r( Ol )...r(ai ) 



a 8 +n 7 o 9 +ra 8 ol0---"5-n7-ns 
x 8 x 9 x 10 ' 

(74) 

with the notation 



r (-j + ni + n 2 + n 3 + n 4 ) r (-n 2 - n 3 - n 4 + n 5 ) r (-j + ni + n 2 + n 6 ) 
^ (ni '-' ni ° )= r(f +m) r(-n 2 -n 3 -n 4 ) r(f +n 1+ n 2 ) 

T + ni + n 5 + n 7 + n 8 ) r (f + ni + n 6 + ra g + m ) ni 
X r(f +ni+n 5 ) r(f + ni + n 6 ) 

The value of G is now obtained as a direct application of GRMT as 

G " rM. .r (o .o)|dc t( A)| '""=-""»n r '-;> < 75) 

where the variables {n*} are solutions of the linear system 





1 























1 


^ 




( »i ^ 




( 




\ 































1 




n* 2 












-i 














-1 








-1 


-1 




n 3 






-a 3 + f 






i 


1 











1 
















nl 






— a 4 












1 































-a 5 















1 






















n% 






-a 6 









-1 


-1 


-1 


1 



















77,7 






-07 +f 
























1 













n% 






-a 8 



























1 










n l 






-a 9 




V 














-1 





-1 


-1 





o J 




V »10 J 




{ 


— aw + -j 


/ 



(76) 



The determinant of this matrix is 1 and the indices n* become 
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ni 



&4 — Q5 — CLq 

, 8 + a 9 + aio, 



0,7 — a 8 — ag — aio, 



no 



3D/2 — a\ — a.2 — &3 
—D + a 5 + a 6 + a 7 + 

-0*5, 

-a%, 

-D/2 + a 8 + a 9 + a 10 , 
-D/2 + ai + a 2 + a 3 , 
-as, 
-ag, 

-3.D/ 2 + a 2 + ft 3 + a 4 + a 5 + a 6 + a 7 + a 8 + a 9 + aio, 
= -a 2 , 



(77) 



that yields the value of the diagram as 



G= (-1)- 



^ ("T — a 89,lo) r (-j — ai 23 ) r — a 2 3456789,lo) T — 0123456789) 

r (ai) r (a 4 ) r (a 7 ) r (a 10 ) r (2£> - ai 23 456789,io) 



r (ai23456789,10 — ^) F (D — a 56 789,lo) T (£> — ai23456) r (-j — CI7) T — a 4 ) 

r (d — a7 89! i ) r (z) — 01234) r — 01234567) r — a 4 5 6 7 89 ,io) 

X ^^—1123456789,10 



(78) 



with the notation 



a-ijk... = &i + cij + dk + ••• 



(79) 



employed above. 

A relevant special case is when all powers of propagators are 1; that is, en = 1 for i = 1, • • • , 10. 
This is given by 



r f _ ir * r(io-^)r(f- 3 )-r(^-9)-rp-6)^r(f-i)- 

U ~ [ L> '3D ' 



r (2D - 10) r (d - tt r - 7) 



(80) 



5 Conclusions 



This paper introduces a technique (GRMT) for the evaluation of a large variety of Feynman dia- 
grams. The advantage over previous methods is that the evaluation of diagrams is reduced to series 
expansions of the integrand, coupled with the solution of a linear system of equations. 

The method is illustrated here in diagrams where the number of series appearing in the process 
is the same as the dimension of the integrals involved. Future publications will describe examples 
where this condition is not present. 
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